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We propose a method for determining ultra-violet divergences in the vacuum energy for systems 
whose spectrum of perturbations is defined through a non-linear spectrum problem, i.e, when the 
fluctuation operator itself depends on the frequency. The method is applied to the plasma shell 
model, which describes some properties of the interaction of electromagnetic field with fullerens. We 
formulate a scalar model, which simplifies the matrix structure, but keeps the frequency dependence 
of the plasma shell, and calculate the ultra-violet divergences in the case when the plasma sheet is 
slightly curved. The divergent terms are expressed in terms of surface integrals of corresponding 
invariants. 

I. INTRODUCTION 

Recent advances in experimental techniques and in the nano-sciences have caused an increase in the interest to the 
theory of the Casimir effect [H, 0|- The quantum field theory (QFT) methods proved to be very efficient in the study 
of the Casimir effect. However, real materials behave rather different from what we have got used to in conventional 
QFT. The properties of these materials change drastically at short distances and/or high frequencies. Numerically, 
the details of this behavior may give a dominant contribution to the Casimir energy [3[ ■ The main focus of the present 
paper will be on the plasma sheet (or plasma shell) model introduced a few years ago by Barton [1, [f| [|| to describe 
certain features of interaction of giant carbon molecules (building blocks of fullerens) with the electromagnetic field. 
The matching conditions for the electromagnetic field become in this model frequency-dependent, and this changes 
drastically the spectrum of the fluctuations. 

The first problem one faces when applying the QFT methods to the Casimir energy calculations is the evaluation of 
divergent terms and renormalization. Though calculations for selected simple configurations may help to understand 
the behavior of the finite part of the Casimir energy, the requirement of selfconsistency (renormalizability) demands 
precise form of the counterterms on considerably more general backgrounds. This can be obtained relatively easy if 
the interaction with the background does no depend on the frequency, but becomes a tough problem if it does. 

To understand the problem better let us start with a frequency-independent case. Consider a system where the 
spectrum of the fluctuations is defined by a linear equation 

Difiu = uj 2 ^ , (f) 

where lu is the frequency of the fluctuations, and I? is a Laplace type operator depending on background fields or 
other external conditions, as, e.g., boundaries. Let p(oj 2 ) be the spectral density of the eigenvalue problem {T]). Then 
one can define the vacuum energy as 

e = -J d{u 2 ){^p{u 2 ). (2) 

(We use the units such that H = c = 1). \i is the lower bound of the spectrum. The eigenvalues lu 2 are supposed 
to be real. The integral ^ is typically divergent at the upper limit and has to be regularized and renormalized. In 
many regularization schemes the divergences are defined by the heat kernel coefficients ak(D) with k < n + 1, with n 
being the number of the spatial dimensions. The heat kernel coefficients are the coefficients in a small t asymptotic 
expansion of the heat trace of the operator D, 

K{D, t) = Tr (e- tD ) ~ t {k - n)/2 a k (D) . (3) 

The coefficients a k are local invariants, which ensures locality of the counterterms, and are known up to a rather high 
number k for a very large class of spectral problems 0, B 0j- Note, that the dependence of D on the background fields 
may be practically arbitrary, as well as the shape of the boundaries. The heat kernel expansion defines therefore the 
ultra violet divergences at the one-loop order for most of the relevant physical problems. 

The situation becomes more complicated if the operator D itself depends on the frequency, i.e. if we are dealing with 
a non-linear spectral problem. If the dependence of D on lu is rather mild, as happens on stationary but non-static 
backgrounds [HI HH, > the divergences can still be expressed through the asymptotics of the heat kernel for D 
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where the dependence on u> is frozen. In time-space noncommutative theories, which also lead to non-linear spectral 
problems, the divergences of the vacuum energy are also expressible through the polynomials of the background 
field, but the structure of these polynomials becomes more complicated [HI, EH- However, it is unclear whether this 
property holds in a more general case. For example, one can calculate the divergences of the vacuum energy in a 
dielectric with a frequenc y-d ependent permittivity, but only for a dielectric body having the form of a ball, or of a 
cylinder, or of a cube flEl. liH fl7j] . 

Can one nevertheless find the asymptotic of spectral density and the ultra-violet divergences for a fairly general 
class of geometries and background fields? The answer is positive, in principle (i.e., modulo technical difficulties). The 
method we propose consists in the following. First one has to freeze the dependence of the operator D on frequency 
and construct the corresponding heat kernel, which may have the form of an expansion in small variations ua of the 
background fields, form of the boundary or boundary conditions, but has to be exact for all values of t and for all 
values of the parameters in D which will later depend on the frequency. From this heat kernel one calculates the 
spectral density for the frequency-independent operator, which, in turn, defines the spectral density for the original 
problem. The high-frequency behavior of the latter defines the ultra-violet divergences of the vacuum energy. 

An important remark is in order. All expansions we shall actually consider are asymptotic expansions in local 
invariants constructed from ua and their derivatives ordered according to the canonical mass dimension. In the 
framework of such an expansion it does not make much sense to compare actual numerical values of the invariants or 
too look for non-local terms. The usefulness of such an expansion is based on the observation that the counterterms in 
QFT arc local. On the other hand, application of the same expansion to calculations of the finite part of the vacuum 
energy is risky, even though the expansions for the spectral densities may formally be integrated as in @. 

We test our method at the example of a slightly curved plasma sheet (more precisely, at a scalar model which keeps 
essential features of the TM model of the electromagnetic field interacting with the plasma sheet). At the end, we 
find the part of the spectral density which has the form of an integral of local quantities quadratic in the extrinsic 
curvature, and which generates logarithmic divergences of the vacuum energy. In the conventional QFT approach this 
divergent part defines local counterterms which must be added to the action. In a more phenomenological approach 
it defines the behavior of the model near the cut off. 

This paper is organized as follows. In the next section we describe the method and study a simple example of 
a frequency-dependent mass term. In section IIIII we discuss the plasma shell model, including the corresponding 
matching conditions, the heat kernel for a frequency-independent auxiliary problem, and the spectral densities with 
their ultra-violet behavior. The last section contains a brief discussion of the results and some future prospects. 

II. GENERAL STRATEGY 

Let us consider a non-linear spectral problem of the type 

D(Q(w),u A )<p = u 2 cp, (4) 

where D is a second order partial differential operator which depends on the frequency u> through a parameter Q(lu). 
The parameter Q may be present in the symbol of D (as, e.g., a mass term), or enter the spectral problem through 
boundary or matching conditions (cf. sec. IIII|) . The parameters ua are supposed to be small. A good example for ua 
is a background field. We suppose that the spectrum of w 2 is bounded from the below for reasonably small ua- Next, 
we assume that for all relevant values of Q there exists the heat trace 

K(D(Q,u A ),t) = TT(e- tD ^^). (5) 

Let us expand K(D(Q,u A ),t) in asymptotic series, 

K(D(Q,u A ),t) = K {0) (D(Q,u A ),t) + K {1) (D(Q,u A ),t) + K {2) (D(Q,u A ),t) + ... (6) 

We shall consider exclusively the case when ua are fields having a positive mass dimension. The expansion ([6]) can 
then be arranged according to the total mass dimension of the fields entering each term, counting also the derivatives. 
I.e., K/q-s does not depend on ua, is linear in u a , K(2) contains the terms quadratic in ua and linear in the 
derivatives of ua, etc. The structure of the expansion ([5]) is similar to the standard heat kernel expansion ([3]), but 
the numerical coefficients in front of local invariants constructed from ua become functions of Q and t. The existence 
of the asymptotic expansion ([5]), and the actual construction of such am expansion, are the most non-trivial and 
demanding parts of our procedure. 

The heat trace can be represented through a Laplace transform of the spectral density ^(A, Q; ua), 



p oo 

K(D(Q,u A ),t)= dXe- xt p(X,Q;u A ), 



(7) 
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where p is the lower bound of the spectrum. The expansion ((6|) defines through the inverse Laplace transformation 
an expansion of the spectral density 

p(X,Q,u A ) = P( )(\,Q,ua) + P(i)(X,Q,u A ) + p(2){X,Q,u A ) + ■■■ (8) 

Next, we have to construct a spectral density corresponding to the initial eigenvalue problem Consider the 
eigenvalue problem 

D(Q,u A )(fQ,x = Xlp q> \ ( 9 ) 

with a fixed Q. The density p(Q, X;u A ) corresponds to this eigenvalue problem, i.e., p{Q, X;u A )dX is the number of 
solutions of ([5]) with the eigenvalues between A and A + dX. Of course, the functions <Pq(w) do solve the equation 
(|4|), but the density is not simply p{Q(tS), uj 2 ; u A ) (with respect to duj 2 ) since the spacing between the eigenvalues 
changes. However, the number of the eigenvalues up to certain to 2 — X remains the same. Namely, if we define a 
counting function 



N(X,Q;u a )= p(u,Q;u A )du (10) 
Jfj, 

then 

N{lo 2 ;u a ) = N(uj 2 ,Q(uj),u a ) (11) 
is the counting function for the non- linear spectral problem (|4]), and 

p{uo 2 ;u A ) = ^Ln(lu 2 ;u a ) (12) 

is the corresponding spectral density (with respect to dui 2 ). This result also follows from a rather complete analysis 
of non- linear spectral problems presented in Refs. [ToL ITU 

The expansion J5J) defines an asymptotic expansion of the spectral density of the non-linear spectral problem 

p(uj 2 ;u A ) = p( a ){uj 2 ,u A ) + p {1) {u) 2 ]u A ) + p( 2 )(uj 2 ,u A ) + ... (13) 

Consider now the limit lo 2 — > oo. According to ([2]) the part of the spectral density which decays as a; -3 or slower 
generates ultra-violet divergences of the vacuum energy. Corresponding counterterms can be isolated in each order of 
the expansion in u Al and the renormalization procedure can be carried out. As in the normal case, the counterterms 
are expected to be local invariants constructed from u A . 

Note, that pu-\ are only parts of the full spectral density. Therefore, p/^ must not be positive. 

Let us give here a short summary of the proposed method. If one likes to evaluate divergences of the vacuum energy 
for a system those fluctuations are described by the equation up to a certain order in u A , one has to make the 
following steps. 

• Consider an auxiliary problem with the dependence of Q on lu frozen and calculate the heat kernel for D up to 
the desired order in u A but exactly in Q and t. 

• Calculate the corresponding spectral density and the counting function N(X, Q;u A ). 

• The function N(uj 2 ,Q(lu);u a ) is then the counting function for our initial problem, and the derivative of 
N(uj 2 ,Q(uj):u a ) with respect to uj 2 gives the corresponding spectral density. The behavior of this spectral 
density for large uj defines ultra-violet divergences of the vacuum energy. 

Note, that the first two steps above are most difficult from the technical point of view, and these steps should be done 
for the auxiliary problem. Therefore, what really counts for the complexity of the method is how D depends on Q 
rather than how Q depends on lu. 

If we are considering only local terms (in u A ) in the expansions above, we would miss non-local divergences it they 
do exist. Non-local divergences arc, however, to exotic to consider this is a serious drawback of the method. 
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A. Frequency-dependent mass term 

As a toy model, let us consider a spectral problem in n — 3 spatial dimensions with the operator D of the form 

D = -df + m 2 + V(x), (14) 

where V(x) is a smooth potential which falls of at the infinity of R 3 , m 2 is a mass term, which will be made frequency- 
dependent later. We shall be interested in the terms depending linearly on the potential. The corresponding part of 
the heat kernel can easily be found (see, e.g., @) 

K {1) ( m 2 ,V,t) = -^^ J d 3 xV{x)e~ mH . (15) 

By using ([7]) and (|10p one immediately finds the corresponding spectral density and the counting function 

P(1) (A,m 2 ;^) = --i I (A-m 2 )- 1 / 2 0(A-m 2 ) / d 3 xV(x), (16) 



8tt 2 

N {1) {X,m 2 ;V) = --^(X-m 2 )^ 2 e(X-m 2 ) f d 3 xV(x) , (17) 



4tt 2 

where 9 is the step function which simply tells that the lower bound of the integration in ([7]) and (|10p is [i = m 2 . Let 
us now consider examples of the dependence of m 2 on u>. 

n m 2 does not depend on the frequency. By expanding (A — m 2 )^ 1 / 2 in (| 16(1 for large A = uj 2 , (A — m 2 ) -1 / 2 ~ 
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1 + ^m 2 uj 3 + . . . we find two divergent contributions to the vacuum energy. The term J d 3 xV diverges linearly, 
and the term with J d 3 xm 2 V diverges logarithmically. This is a well-known result. 

b. m 2 = aoj 2 . Taking A = uj 2 we obtain (A — m?) = (1 — a)ui 2 . Therefore, for a < 1 all spectral functions are as 
in the case above with m 2 = up to a multiplier of y/1 — a. For a > 1 the spectral density is non-zero for negative 
uj 2 meaning that the system becomes unstable. The same conclusion may be reached by considering the eigenvalue 
problem uj 2 ip u = (aw 2 — d 2 + V)<p u . 

c. m 2 = juj . For 7 > the spectral density is non-zero for < u> 2 < 1/7, i.e. there are no high frequencies 
in the spectrum at all. Consequently, there are no ultraviolet divergences proportional to the potential. For 7 < 0, 
instabilities appear due to the imaginary frequencies with ui 2 < — 1/|7|. 



III. CURVED PLASMA SHEET 



A. Matching conditions 

The plasma shell (or plasma sheet in the noncompact case) model was introduced by Barton in Refs. 4, 5] where also 
the Casimir energies were explored. The Casimir force between the plasma sheet and other surfaces or molecules was 
calculated in [1 81 ] by using the Lifshitz formulae. Further spectral properties of the plasma shell model were studied 
for flat [l9| and spherical [2(| geometries. In both cases, the TM mode of the electromagnetic field causes more 
difficulties than the TE mode. The TE problem is isomorphic to the scattering on a repulsive frequency-independent 
delta-potential, which is a rather well understood problem. The study of the heat kernel expansion and thus of the 
divergences of the vacuum energy for the delta-potential on various surfaces started in Ref. [21( with the case of a 
sphere. Then the heat kernel for an arbitrary shape of the surface in an arbitrarily curved space, with generic bulk 
potential and gauge fields was obtained [12] • The spectral problem for TM mode is known to be not elliptic [lij H(| • 
Thus even the definition of basic spectral functions become problematic. 

Here we consider exclusively the TM mode, which, for a spherical shell placed at r = tq may be reduced to a scalar 
field ip with the following matching conditions [J] 

discont {npu) = -■^■(9 r r^), 

discont (d r r<fu) = , (18) 

where discont (/) = f(r = ro + 0) — f(r = tq — 0). q is a constant depending on the properties of the material. Outside 
the shell the modes ip u satisfy the Helmholtz equation 



(V 2 +w 2 )^ =0. 



(19) 
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In the limit ro — > oo the matching conditions (JTHJ) reproduce that of the fiat plasma sheet 0, . 

To prepare for the discussion of the next subsection we need a bit more elaborate notations. Denote by M + the 
exterior of the shell, and by M~ - the interior of the shell. Let E be the (spherical) interface of M + and M~ , and let 
tv 11 be a unit normal to E pointing inward of M^. Denote corresponding covariant derivatives by V^. When acting 
on scalars, Vlf = ct, and V~ = —d r . Then we can rewrite (fT8l) as 



V n <p- ) \ -Q Q-tr J \ tp- 



(20) 



where 



QH = u 2 /(2q), a = l/r . (21) 



The dependence of ip on lu is suppressed. 

We shall consider the case when the surface E has the topology of a plane but is slightly curved. M + and M~ are 
two parts of the bulk manifold separated by E. The non- flatness of E is measured by the extrinsic curvature tensor 
L a b which has two indices corresponding to the directions tangential to E. Obviously, the extrinsic curvature of E 
considered as a part of M + is minus the extrinsic curvature of E as a part of M~ , L^ b = —L~ b . We shall consider 
a scalar field with the matching conditions l|20p on E, where, according to (|21[) . Q is independent of the extrinsic 
curvature, while a is supposed to be of the order of \L a b\, but independent of the frequency. Therefore, Q will be 
treated non-perturbatively, while only the leading orders of a and L ao will be retained in all spectral functions. 

Note, that for an arbitrary curved boundary the electromagnetic field cannot be separated in the TE and TM parts. 
Therefore, the model we are going to analyze is an approximation to real physics. We hope that the most essential 
features of the problem are nevertheless preserved by this simplification. 



B. Heat kernel for the auxiliary model 

The auxiliary model is obtained by simply freezing the frequency dependence in the matching conditions (|20j) , i.e. 
by assuming that Q is some constant. We are going to evaluate the localized heat kernel 

K(J,D,t) = Tt(fexp(-tD)) (22) 

where / is some smooth localizing (or smearing) function, and D is the standard Laplacian, D = — V 2 . 

The method we are going to use is (an extension of) the method of Gilkey [23|, [24[ , see @, H, Q for an extensive 
overview with further explanations. The method is based on a somewhat paradoxical observation that solving a more 
general problem may be much easier. We introduce an arbitrary curved metric gij on M and generalize D to be an 
arbitrary scalar Laplacian 

D = -(g»V i V j + E), (23) 

where E(x) is a potential. We suppose that the metric is smooth across E, so that L^ b = —L~ b . At some stage, when 
discussing the doubling trick, see Eq. (|28p below, we shall need more general matching conditions 



v n¥ >+ \_ s++ s+- \ /V 
v n <p- I { s-+ s— ) { <p- 



(24) 



depending on four functions S f±± . The functions 5' ±± have a positive mass dimensions and are, therefore, natural 
parameters of a perturbative expansion. However, S ±± — does not correspond to a free propagation, but rather 
to two disjoint regions with Neumann boundary conditions on the surface which separates them. On a side note, 
we remark that this explains the failure of the multiple reflection expansion in the case of these matching conditions 
[25l | . Throughout this subsection we shall suppose that the dimension of M is arbitrary, n = dimM. The last (nth) 
coordinate is supposed to be normal to E, so that the notations introduced above for the normal and tangential 
coordinates remain consistent. The spectral geometry of Laplace type operators with such matching conditions has 
been analyzed in (26j . Basing on the analysis of that paper we can write the following expression for the heat kernel for 
the conditions (|20|) which is non-perturbative in Q and contains other invariants up to the canonical mass dimension 
two: 

D > *) " ( 4 ^)(n-l)/2 £ ^'^ ^ faWf + - L * b )( L t b L ao)f + foWL+a L~ a ff 

+0s(z)(f+ - /")(£+ - L- a ) + fc{z)\{f+ n + /"J + p 5 (zHL+ a - L- a )f + fc(z)v(f+ - /") 
+l3 7 (z)Ef + s (z)Rf + /3 9 (z)R anan f + l3 10 (z)a 2 f)] + K vol (f, D, t). (25) 
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Here 0i{z) are unknown function depending on Q through the dimensionless combination 

z = 2Qt 1/2 . (26) 

(The results of the calculations are summarized at the end of this subsection). Q is supposed to be constant, but a 
is allowed to depend on the coordinates on S. The volume part cannot depend on Q or a and is thus given by the 
standard expression 

K vol (f, D, t) = J d n x^gf(l + t(E + iZ/6)). (27) 

We stress again that only the invariants of canonical mass dimension two or less are being kept. Only a few functions 
Pi will be needed to analyze the curved plasma sheet model. Others are included for technical reasons. 

Some explanations regarding the formula (|25p are in order. Rijki is the Riemann tensor for the metric gij , and R is 
the corresponding scalar curvature, h is the determinant of the metric induced on £. All tensor indices are taken in 
a local orthonormal frame, so that there is no distinction between upper and lower indices. n ± is a unit vector along 
a normal geodesies pointing inward M . The indices a, b correspond to the vectors tangential to S. Summation over 
the repeated indices is understood. Semicolon denotes the covariant derivative, e.g., /.+ n = V+V+/ + . The smearing 
function / and the metric are smooth across E, consequently, /.+ = — /.", L+ fe = —L~ b , /.+ n = f~ n . Therefore, the 
notations used in (j25|) are redundant, but they facilitate the comparison to Ref. [26j and allow to keep track of the 
symmetries more easily. Eq. (|25|) contains all possible invariants up to the mass dimension two. Note, that the heat 
kernel must be invariant under interchanging the roles of M + and M~ . Since a is antisymmetric under M + «-> M~, 
this excludes the surface integral of a as well as all other (non-vanishing) invariants of the mass dimension one. Since 
a is defined on the interface only, one is not allowed to differentiate it with respect to the normal vector. Finally, 
Rabab is not an independent invariant. 

Here we like to reiterate the remark made already a couple of times above. The expression (|25|) is an asymptotic 
expansion in local invariant ordered according to the canonical mass dimension. It has nothing to do with actual 
(numerical) smallness of corresponding terms. 

Let us start the calculations of the unknown functions Pi{z). By standard arguments (see, e.g., section 4.1 of Ref. 
Q) one can show that (3i{z) do not depend on n, i.e. the whole dependence on the dimensionality of the manifold 
resides in the pre-factor (^Tri) - ^ -1 " 2 . This fact is extremely important for the methods we use here. As we shall 
see below, the so called conformal relations yield some relations on 0i with coefficients which are linear in n. Since pi 
themselves do not depend on n, each of the conformal relations in fact produces two independent conditions on /Vs. 

The functions (3q, P4, (3§ and P\q can be computed by using Lemma 4.1 from [26j which relates the heat trace on 
a manifold with boundary with Robin boundary conditions to the heat trace on a "doubled" manifold with transfer 
conditions on the interface. This Lemma states the following. Let Af ± = M° be a Riemannian manifold, and let 
Z) ± = D° be a scalar operator of Laplace type. Fix an angle < 9 < n/2. Let S ++ and be arbitrary. Set 

S-+ :=£+-, 

S— := S ++ + (tan - cot 0)S + ~ , (28) 
SV := S ++ + tnn6S + -, 
:= S ++ - cot 6S+-. 

Then 

K(f, D, t) = K(cos 2 0f+ + sin 2 Of-, D% t) + K(sm 2 0f+ + cos 2 9f~,D%t). (29) 

On the left hand side of the equation above we have the heat kernel for the operator D on M subject to transfer 
conditions on the interface, while on the right hand side there are two heat kernels on for the operator D° with 
Robin boundary conditions 

(V„</> + S^ldMO = and (V„-0 + S^)ip\g M a — , (30) 

respectively, is a restriction of / to M^. 

To control the invariants appearing with /3q, /34, Pe and /3io it is enough to consider the case of constant a. Then, 



(7 



t a n0=^ + J- 2+ l = l + - + w+ O(**) (31) 
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and 

2Q ' ~ v ~ " ~ v ' 2Q 



^ = -^ + 0(ct 3 ), 5v, = 2Q+^- + 0(a 3 ). (32) 



Another ingredient we need, is the heat kernel for the free Laplacian A = — df subject to the Robin boundary 
conditions at a plane boundary with a constant S which was derived in [27j 

K(f,A s ,t)= - J 1W9 / <r- 1 xY(Yfto(st 1 '*y+ k .TT 1 , +/fr>f/ a XT* , Y (33) 

where is the pth normal derivative of /. We do not write here a trivial volume part of the heat kernel. For each 
p one can also write a closed expression in terms of the error function 

erf (z) = 1 - erfc (z) = -= / e^du. (34) 



/7T 



For example, for p = we have 

A - " - mh* L [eS " ert (s ' ,,2) + eS " - 5' • «»' 

In the formulae above, Eqs. (|33p and (|35p we corrected the ^-independent terms of the corresponding formulae in 
[27| . Such terms, being of the order of "zero reflections" were outside of the focus of the paper (27j . 

Let us stress, that as noted in the expression (|35[) is exact, i.e. it takes into account the possibility of non- 
analytic contributions which may not be seen at power-series expansions, see also [28l |29|. 

By collecting ([28]) - (f35]) together, we obtain, 

A,(z) = i[ e * 2 (erf(z) + 1)] 

2 to r(| + i)' 
1 00 1 

f3 6 {z) = 2fa{z), 
P 10 (z) = 2p (z). 

Some terms in the expansions above can be checked against the corresponding terms in [261 ] - These expansions are 
rapidly convergent for positive z. Therefore, it does not play a role whether one uses closed formulae in terms of the 
error function, or an expanded form. Using expanded forms is a little bit easier. 

To define (3q and fa we use the product Lemma (see, e.g., sec. 4.1 of Ref. H) which states that if M = M\ x M 2} 
and also D = D\ ® 1 + 1 ® D2, with D\ and D2 acting independently on functions over Mi and M2, the resulting 
heat kernel factorizes into a product of the heat kernel for D\ and the heat kernel for This means that if 

f(x\,X2) = fi{xi)f2{x2), where x\ and x 2 are coordinates on M\ and M 2 , respectively, 

K(J, D, t) - K(fi, D u t)K(f 2 ,D 2 , t) . (37) 

Let us take M 2 being a curved manifold without a singular surface, and M\ being a flat manifold with a singular 
surface and define some transfer conditions on it. Let D\ and D 2 be scalar Lapalcians, and let E be independent of 
the coordinates on M\. Then, to the linear order of E and the scalar curvature R 

K(f 2 ,D 2 ,t)= (47rf) *_ na/2 / dx 2 ^T 2 {E + \R). (38) 

Then, by comparing corresponding terms on both sides of (|37p . we obtain 

/3r(«) =/%(«), fo(z) = ±Po(z). (39) 
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The remaining functions [3\, 02, Ps, /3s and /3g will be computed by using conformal properties of the operator D. 
One can show that 



d_ 

de 



e=0 



Tr [exp (-e- 2ef Dt)] = -M—Tr [/exp(-tD)]. 



(40) 



It is easy to figure out (or to look up in [7, 8, 9]) how the geometric quantities entering the operator D, Eq. (f23j) . 
transform the local scale transformation D — > e~ 2e f D. In particular, the transformation of the metric is the standard 
Weyl rescaling, — > e 2(L ^ gij- The transformation of other quantities, as the connection and the potential, are less 
standard, as they contain specific contributions which ensure the homogeneous transformation law for D. For our 
purposes, it is enough to consider the transformations vanishing on S, /|s = 0. Under this assumption, the variations 
of relevant geometric invariants on the interface E read 



d 
1 
d 
1 
d_ 
1c 
d_ 
de 



E = \{n 2)(/+„ + /"J - \{n - 2)(Z+ - L" )(/+ - /") , 



R=-(n- l)(/+„ + /.;„) + -(L+ - L" )(/+ - /-) 



R, 



"ii".n— 2 ( n l)(/;nn + /;mi) ^ (^oa ^ao)(/;n f;n) 



(41) 



c=0 



We must keep the matching conditions invariant under the local conformal (scale) transformations. This is achieved 
by transforming a to compensate the variation of the connection in V n , 



d_ 

de 



e=0 



<r=\(n-2)(f+-f-) 



(42) 



(cf [2a]). Q does not transform. 

Let us collect the terms with | (/."£„ + /.^„) in (|40|) . This gives the equation 



1 



(n - 2)/5 7 (z) - 2(n - l)/3 8 (z) + (n - l)/? 9 (z) = [(n - 3) - zd z ]^{z) , 



which, with the help of ([59^1 . may be simplified to 

1 



G 



(n - 4)/3b(«) + (n - !)#,(*) - [(n - 3) - z&J&O*) . 



The functions Pi(z) do not depend on n. By taking n = 1 we obtain the identity 
which may serve as a consistency check. Taking n — 4 yields the value of /3g, 

Next, lets us collect the terms with cr(/.+ — /.~) in (|40p. This yields the equation 



(1 - n)f3 5 (z) + \{n- 2)/3 10 (z) = [(n - 3) - z^]/3 6 (z) . 
Again, n = 1 provides a consistency check. For n — 2 we obtain 

AW = [i + zdMz) = \ E^^tj^j • 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 
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The only structure which has not been used yet in the conformal relations is (/.+ — f. n )(L+ a — L aa ). It gives the 
condition 

i(2 - n)(3 7 {z) + i(n - l)fo(z) ~ J/3 9 (z) + \{n- 2)/3 5 (z) + 2(1 - n)(i 2 {z) - 20i(z) = [(n - 3) - zd 9 ]fa(z) (49) 



or 



^(4 - n)p (z) - l -(3 9 (z) + l -(n - 2)fc(z) + 2(1 - n)p 2 (z) - 2fr(z) = [(n - 3) - zd z ]p 3 (z) 



(50) 



Since the functions fli(z) do not depend on n, this equation is equivalent to two independent conditions. However, 
we still have three unknown functions f3\, (3%, and /?3. 

To get one more relation between the unknown functions we shall use the properties of the heat kernel under local 
scale transformations when the smearing function is also transformed. It is a purely formal computation [231 . [24| to 
show that 



d_ 



l - 2 (e- 2 <fF,e- 2 ^D) = Q, 



(51) 



-=o 



where, similarly to ((3]), ak(F, D) are the coefficients in an expansion of the smeared heat kernel K(F, D, t). One has to 
note that for each dimension n of the manifold the equation (|5T|) restricts just one heat kernel coefficient, in contrast 
to (|40p. Let /?* be the coefficient in front of z k in the Taylor expansion of fti(z), 



(3 i (z) = J2z k $. 



(52) 



fc=0 



Since z ~ t 1 ^ 2 , cf (|26| . each (5% contributes to I.e., Eq. (f5Tjl gives restrictions on /3 fc in the dimension n = k + 5. 

Lets us collect the terms in ([51"]) proportional to (/.+ — f^)(F+ — F.~). Since now the smearing function is also 
transformed, such terms may also appear from conformal transformations of the second normal derivative of F: 



d_ 

de 



(F^ 



F~ ) = --(f + - D(F + - F~) + 



(53) 



(cf. Ref. [2J|). Dots denote other terms, e.g. ~ f-,, m F, which we do not use in this calculation. We have, 

0. (54) 
Next, we remember that n = k + 5 and use (f3"6")) to obtain 



^-(n-l)/^ + i(n-2)/3 6 fe 



64 T(|+3) ' 

Now we are ready to compute the remaining functions. Taking n — 1 and n = 2 in Eq. (|50[) gives (3\ and /?2, 
respectively, 



1 OC 

-T 



192 



fe=0 



-. oc 

A2V ' 384 ^ 



fc=0 



fc (fc + l)(2fc + l) 

r(f + 3) ' 

fc (fc + l)(5fe + 13) 

r(f + 3) 



(56) 



For the convenience of the reader we summarize the results of calculations of this subsection in Table U by providing 
the references to corresponding equations. One can easily check that that all 0i(z) are consistent with the results of 
|2(| where several leading terms of the z-expansions for each (3 were calculated. 

The last step is to sum up the z-expansions for selected functions f3i(z) which are going to use below in sec. IIIIDl 



4 dz 



-(e z (erf(z) + l)-l) 



Hz) 



l 

12 
5 

48 



192^' 
3^' 



(57) 

(58) 
(59) 
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TABLE I: Summary of the calculations of /3i(z). 



Pi-. 


/3o, At, Pe, Pio 


Pr, Ps 


p 9 


Ps 


Ps. 


Pi, P2 


Equation: 






gel) 









where 



0(z) = v- 
az 



1 



2z 



-r e 2 (erf (z) + 1) - 1 - -= - z 



(60) 



Note, that in principal the power series expansion, even though they are convergent for all positive z, define the 
functions up to non-analytic terms which are exponentially small at z — > 0. Nevertheless, since the expression for the 
Robin heat kernel (|35[) is exact and does not include non-analytic parts, the functions /Jj with i = 4, . . . , 10 also do 
not include non-analytic parts, and their expressions in terms of the error function are also exact (though not all of 
them are presented explicitly here). However, to obtain f3$ the use of an expanded form was essential to write down 
the conformal relation (|54[) . Therefore, (3^{z) in principle can contain a non-analytic part. Let us assume that the 
non-analytic part of (3z{z) is zero, which is a very natural assumption. Then, there is no non-analyticity in (3\ and (3% 
which are defined by the other functions through algebraic equations, and f3i and 02 are restored in a unique way as 
written above. 



C. Spectral densities: plane case 



We have completed the calculation of the heat kernel, and it remains to calculate the spectral densities. However, 
there is a subtlety in application of our methods in the present case, which is related to the treatment of unstable 
modes with negative u> 2 . To understand the problem we start with the case of a plane plasma sheet. 

One can easily see that flat (a = 0, L^ b = 0) auxiliary problem (with Q positive and independent of us) admits 
surface modes which decay as e -2 ^"' for positive x n and as — e -2 ^" for negative x n . The Helmholtz equation 
(1191) gives the dispersion equation 

oj 2 =P 2 -(2Q)\ (61) 

where p is the momentum parallel to the interface S. This dispersion relation corresponds to the excitations localized 
at the surface of the interface and having a negative mass squared —AQ 2 . Due to the presence of such tachionic 
excitations, the heat kernel has a part which grows for large positive t. This part is defined by the growing part of 
(3o(z), cf. Eqs. (|2"5)) and ([55]) . In n = 3 dimensions it reads 

K tachioa (f, D, t) = ^-e 4Q2t [ d 2 xf . (62) 

For / = 1 the heat kernel above is divergent, which is a consequence of the translation invariance of the problem. 
Due to this invariance, the relevant quantities are densities per unit area of the interface, as, e.g., the surface Casimir 
energy density. Therefore, we omit the integral of the smearing function in Eq. ()62|1 thus obtaining a surface density 
of the heat trace. The corresponding spectral density and the counting function 

Ptachio„(A, Q) = j-#(A + 4Q 2 ), (63) 
^tachio„(A,Q) = -^(A + 4Q 2 )0(A + 4Q 2 ), (64) 

47T 

also have the meaning of surface densities. By using Eqs. (fTTT) . (|T2]) and (|21"j) . one obtains the part of the spectral 
density for the nonlinear spectral problem, which is generated by the tachionic part of the auxiliary problem: 

Ptachion(w 2 ) = [l + 9 (iO 2 + \ . (65) 



This spectral density has two branches. One corresponds to real frequencies, ui 2 > 0, while the other - to imaginary 
frequencies, lo 2 < —q 2 . 
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The flat plasma sheet model can also be analyzed directly. There is a mode, the so-called surface plasmon, which 
falls exponentially away from £ with the cigcnfrcqucncics 19] 

<. = § (V¥+Ap~ 2 -q) ■ (66) 

The spectral density corresponding to the surface plasmon is simply the density of the momenta p parallel to £, i.e., 
(2ir)~ 2 d 2 p, which after integration over the directions of p becomes (4:Tr)~ 1 d\p\ 2 . To obtain the spectral density in 
terms of to 2 one has to change variables in the previous expression which yields 

After some simple algebra the density (I67p can be shown to reproduce the positive branch of the spectral density 
(|65|) . One can also demonstrate that the negative branch of (|65|) . u> 2 < —q 2 , correspond to the solutions which grow 
exponentially away from E and thus do not belong to the physical spectrum. In our approach, the negative branch 
must be excluded by hand. Obviously, in one perturbs the plane plasma sheet by making is slightly curved, this 
cannot result in the eigenvalues far away from the spectrum of non-perturbed problem. 



D. Spectral densities: curved case 



Let us now proceed with the case of a curved plasma sheet. Let us restrict ourselves to the case of flat ambient 
space with no external potentials, so that Rijki = E = 0. Next, we put the smearing function / = 1. If we suppose 
that the extrinsic curvature and a are well localized, the surface integral in (125|) remains convergent (otherwise, we 
can consider the surface density of the vacuum energy as described in the subsection above) . Let us remind, that a is 
supposed to be of the same order as the extrinsic curvature. Therefore, we arrive at the following second-order heat 
kernel 

K {2) (D(Q), t) = j^J^ d 2 xVh[p 1 (z)(L+ b - L- b )(L+ b L~ b ) + 2 (z)(L+ a L 

+P 5 (z)a(L+ a - L" ) + P 10 (z)a 2 ]. (68) 

(the 1st order heat kernel Kn\ vanishes identically). The functions which remain relevant for the calculations of this 
subsection are /3%, 02, /3s and P\o- The numbering conventions do not look very natural, but, unfortunately, any 
change of notations would inevitably bring in misprints. The present author kindly asks the reader for understanding. 
Let us start with the contribution of /3io, cf. Eq. ([56"]) . By using the identity 



0(1 J 



5 4 « 2 *erfc {2Qt 1 ' 2 ) = - / dAe" At ^- arctan 
7T J dA 



2Q_ 

A 1 / 2 



(69) 



we arrive at the following expression for the counting function 



N 1Q (\,Q) = ±- I d 2 xa(x) 2 -N 10 (X,Q), 



4tt 



E 



I' ( 2Q 



N W (X, Q) - 29{\ + AQ Z ) + - arctan I ^ j . (70) 
By using Eqs. (|12[) and (|2Tj) one obtains an expression for the spectral density of the non-linear spectral problem 

Pio(w 2 ) = — j d 2 xa(x) 2 ■ p 10 (uj 2 ) , 

p w {u 2 ) = 28{u 2 )+ q (71) 

TTU>{LU Z + q ) 

We dropped another delta- function, S(u> 2 + q 2 ), since it corresponds to a non- physical part of the spectrum in accor- 
dance with the previous subsection. The part of pio(w 2 ) which decays as w~ 3 at large u>, namely 

P^(- 2 ) = ^~J d 2 ^ 2 (72) 
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generates a logarithmic divergence in the vacuum energy p]). 

Let us now proceed with contributions of other invariants appearing in (|68[) . The inverse Laplace transformation of 
the functions /3s, /3x, P2 (see (|57| - ([6010 may be performed by using e.g. the textbook [33]. The resulting expressions 
for the spectral densities are becoming rather lengthy, and thus we keep track of divergent contributions only. They 
read 

/> dh V) = ibiS/ xVl \-k [L « b - L « b){L °» - L «» ] + ~ L « J2 + a[LL L - ] + ff2] ' (73) 

where we have also included (|72[) . This is the main result of this paper in what concerns the curved plasma sheet 
model. 

The formula (|73j) defines all divergences of the vacuum energy in in the plasma sheet model up to the quadratic 
order in the extrinsic curvature. We like to remind that in sec. IIII Al we argued that a has to be considered as 
being of the same order as L a b- We see, that there are no divergences in the linear order, and that all divergences 
are logarithmic. The expansion in the extrinsic curvature is an asymptotic one. Although it is enough to state the 
presence of divergences of certain kind and the necessity of counterterms, it is not enough to use the spectral densities 
for calculations of the finite part of the vacuum energy, as it is not clear how large or small the higher order corrections 
could be. 

In principle, one cannot totally exclude that some other divergences will also appear in the case of curved plasma 
sheet. Here we mean higher orders of the extrinsic curvature, or with derivatives acting on L a b or a. However, we 
find this very unlikely. Higher order invariants are usually less divergent, and the quadratic terms already diverge 
only logarithmically. 

It is instructive to compare our results to the ones obtained in (20l | for the spherical plasma shell. Since the volume 
of a two-sphere of the radius ro is proportional to t-q, and the extrinsic curvature is inverse proportional to ro, the 
expression (|73p does not depend on the radius. The paper [20j ] found linear divergences which do not depend on 
the the radius and, therefore, correspond to the extrinsic curvature squared. (One should take into account that 
our conventions for the heat kernel coefficients differ from the ones used in [20j.) This means that the divergences of 
plasma model are sensitive to the topology of the interface surface £ in contrast to (most of) the frequency-independent 
problems. This is an expected result, since to derive the part of the spectral densities which generates divergences we 
needed the heat kernel for all values of t, and not only the small t asymptotics. 

From the spectral densities we have calculated above one can derive the the heat kernel coefficients for the non- 
linear spectral problem. This step is, however, unnecessary in our approach since the divergences have already been 
calculated. 

Let us remind that the electromagnetic field may be reduced to scalar fields only for some specific geometries 
(including flat and spherical ones). Therefore, the problem studied in this section is an idealization of the full 
problem. However, the main complication from considering the electromagnetic field directly is a rather nontrivial 
matrix structure of the matching conditions, which makes the calculations more involved, but not to the point where 
the technical difficulties become unresolvable. 



IV. CONCLUSIONS 



In this paper we suggested a method to evaluate the ultra-violet behavior of spectral density and resulting diver- 
gences of the vacuum energy for the problems with generic non-linear dependence on the spectral parameter. The 
divergences are represented in a form similar to the heat kernel coefficients, i.e. through integrals of local invariants 
associated with the problem. The method, however, requires considerably more information on the spectral properties 
of the operator involved that just the asymptotic expansion of the heat kernel. One needs the whole heat kernel for 
all values of t. Also, the parameter which depends on the frequency has to be treated non-perturbatively. 

The method was tested at the plasma model of Barton, for which the divergences of the vacuum energy for of a 
slightly curved plasma sheet were calculated up to the second order in the extrinsic curvature. These divergences 
appeared to be logarithmic, and they depend crucially on the topology of the interface surface. This fact has certain 
implications for the plasma sheet model. Clearly it implies that a renormalizable model (in the sense of [20j ) must 
contain some topology-dependent contributions perhaps corresponding to global degrees of freedom. 

We expect that the method can be considerably improved by using recent advances in related areas. For example, to 
calculate the heat kernel one can employ the worldline formalism 1311 . Recent advances in calculation of the Casimir 
force between bodies of various shapes may also be useful (cf. [32l l33l. [HI, HU, [H, H3, HI] and references therein). 
We also like to note some similarities to the method of Ref. (39| applied to one- dimensional systems with arbitrary 
dispersion relations. Finally, we expect that with some modifications the methods proposed above may be a ppli ed to 
such complicated problems as boundary conditions depending on the tangential momentum studied in Ref. [401 ]. 
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